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Abstract 

The quark mass function in QCD is revisited, using a gluon 
propagator in the form + m^) plus 2fj?/{k'^ + m^^, where the 

second (IR) term gives hnear confinement for nig = in the instan- 
taneous limit, fi being another scale. To find we propose a new 
(differential) form of the Dyson-Schwinger Equation (DSE) for S(p), 
based on an infinitesimal subtractive Renormalization via a differen- 
tial operator which lowers the degree of divergence in integration on 
the RHS, by TWO units. This warrants T,(p — k) Ti[p) in the in- 
tegrand since its ^-dependence is no longer sensitive to the principal 
term {jp — k)"^ in the quark propagator. The simplified DSE (which 
incorporates WT identity in the Landau gauge) is satisfied for large 

by E(p) = i;(0)/(l + /3p^), except for Log factors. The limit = 
determines Sq.A third limit = —rn^ defines the dynamical mass 
mo via 5](«mo) = +mo. After two checks (/^r = 93 ± IMeV and 
< qq >= (280 ± 5MeVf), for L5 < /? < 2 with Sq = SOOMeV, the 
T- dependent DSE is used in the real time formalism to determine the 
"critical" index 7=1/3 analytically, with the IR term partly serving 
for the H field. We find Tc = 180 ± 20MeV and check the vanishing 
of U and < > at T^. PACS: 24.85.+p; 12.38.Lg; 12.38.Aw. 
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1 Introduction 



QCD, as the queen of strong interaction theory, hes at the root of a whole 
complex of strong interaction phenomena, ranging from particle physics to 
cosmology Its principal tool is the quark mass function, termed S(p) in the 
following, as a central ingredient for the evaluation of a string of QCD pa- 
rameters whose primary examples are the pion decay constant and the quark 
condensate. The thermal behaviour of the latter in turn has acquired consid- 
erable cosmological relevance in recent years in the context of global exper- 
imentation on heavy ion collisions as a means of accessing the quark-gluon 
plasma (QGP) phase [1- 4]. It is therefore essential to have on hand a re- 
liable E(p) function in a non-perturbative form as a first step towards the 
evaluation of these basic QCD parameters. In this respect, QCD sum rule 
(SR), attuned to finite temperatures [2] have been a leading candidate for 
such studies for a long time, using the FESR duality principle [5], as well as 
a variational approach via the minimum of effective action up to two loops 
(Barducci et al, [la]) to determine the mass function. An alternative ap- 
proach has been the method of chiral perturbation theory [6] with the pion 
as the basic unit in preference to quarks. Now a standard approach to QCD 
is via the RG equation for the /3 function in the lowest order of g which yields 

a,(g2) = 27r/[91n(Q/AQ)] 

with 3 flavours, Aq being the QCD scale parameter [7]. Unfortunately the 
higher order terms in g are not particularly amenable to the simulation of 
non-perturbative effects. On the other hand, the Dyson-Schwinger Equation 
(DSE) which may be regarded as the differential form of the minimum prin- 
ciple of effective action [8] , offers a more promising tool which has often been 
used with the standard o.g.e. in the rainbow approximation [9], but can be 
improved to incorporate gauge invariance so as to satisfy the W-T identity in 
the "dynamical perturbation theory" (which ignores cris cross gluon lines in 
the skeleton diagrams) with little extra effort, as first shown by Pagcls-Stokar 
[10]. In this paper, we shall use the same approach, but explicitly add an 
extra, non-perturbative, term to the one-gluon-exchange (oge) propagator 
for a quicker simulation of the infrared (IR) regime, so that both together 
act as the 'kernel' of the Dyson Schwinger Equation (DSE) [11]. Thus the 
total gluon propagator is given by 
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where is a scale parameter corresponding to the (hadronic) GeV regime, 
(whose value will be left undetermined till later), and rrig is a (small) gluon 
mass a non-zero value for which can be motivated from several angles, a no- 
table one being the 'Schwinger mechanism' [12] as explained in the Jackiw- 
Johnson paper [13]. A second motivation was highlighted by Cornwall et 
al [14], in the context of their approach to a more compact realization of 
gauge invariance via the so-called 'pinch mechanism' [14]. Yet a third mo- 
tivation which is especially relevant in the present context of a temperature 
dependent DSE, comes from according it a 'Debyc mass' status, running 
with the temperature [15]. A non-perturbative gluon propagator (with har- 
monic confinement) was employed in [16-17] as a kernel of a BSE for the gg 
wave function for the calculation of glueball spectra , on similar lines to qq 
spectroscopy [18] Alternative BSE treatments for glueballs also exist in the 
literature [19]. In this paper, the IR part of Eq.(l.l) has a dual role: 1) 
to serve as a more efficient simulation of the non-perturbative effects on the 
mass function E(p) ; and 2) a partial simulation of the external magnetic field 
effect, as an alternative to small non-zero masses of 'current' quarks [20,15]. 
With a non-perturbative sohition of the DSE, we are primarily concerned 
with chiral symmetry restoration at a critical temperature T^. To that end 
we shall be interested in the T behaviour near the critical Point Tc, rather 
than as an expansion in powers of near T = [6]. Note however that 
linear confinement {-^ r) corresponds to rrig = 0, via the second term in (1.1), 
in the (3D) instantaneous limit t = 0, so that deconfinement competes with 
chiral symmetry restoration with a propagator like (1.1). We shall not pursue 
this aspect further, although we note that deconfinement has been claimed 
to occur at a lower temperature [21] than chiral symmetry restoration. 



1.1 Object and Scope of The Paper 

The central object of this paper is a determination of the mass function 
E(p) non-perturbatively in the intermediate momentum regime with the help 
of the gluon propagator (1.1) that covers the IR regime. This is sought 
to be achieved via a (new) differential formulation of the DSE based on a 
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subtractive form of renormalization that is particularly convenient for a DSE 
type equation. A second object is to apply the T,{p) so determined, to two 
basic quantities, < qq > and fT^jS^nd express them in an analytic form, so 
that their T- dependent generalizations may be achieved analytically too. 
A third object is to generalize the DSE to a T-dependcnt form, so as to 
obtain an equation for the T-dependent mass function rrit, with a focus on 
its critical index 7 associated with the critical temperature Tc, so as to gauge 
the role of the IR term vis-a-vis small current masses to simulate the H- 
field effect [20, 15]. Further, while in the conventional methods [20, 15], the 
various thermodynamic quantities are derived from a central quantity like 
the free energy [15], or equivalently the effective potential [20], and taking 
appropriate derivatives, the plan adopted here is to focus on the DSE itself as 
the principal form of dynamics, with rrit as a natural order parameter. Due 
to the unconventional nature of this approach, this part of the exercise is 
still preliminary, with only one critical index identifiable with rrit determined 
so far, while other indices [20] are left for later studies, within the DSE 
framework. 

In Sect. 2, we formulate the DSE for E(p) in an (infinitesimal) form of 
(subtractive) renormalization which yields a non-linear second order differ- 
ential equation for this quantity. The dynamical mass rriQ is defined as the 
pole of Sp{p) at ij.p = —rriQ, and hence corresponds to the solution of the 
equation for S(imo) = mo- Although in principle a mass renormalization 
factor comes according to rules [7], the condition E(imo) = mo ensures 
that this factor is effectively unity, provided the dynamical mass is employed 
for the propagator at its pole. As for the quantity S(0), we shall designate it 
as the constituent mass. For the solution of the resulting DSE, three crucial 
check-points are = 00; = 0; = —m^ which control the structure 
of S(p). The simplest ansatz consistent with a - like behaviour in the 

= 00 limit, as demanded by QCD, is Eq/ (1 + l^p"^) [1], the only precaution 
needed for a consistent solution being a constant with its argument fixed 
in advance at a certain specified value. This form has good analytical prop- 
erties for large space-like momenta, but it implies that the dynamical mass 
E(iTOo) exceeds the constituent mass E(0). 

For a basic test of this structure, we choose in Sect. 3, two key items i) qq 
and ii) whose derivations are sketched in Appendices A and B respectively 
in an analytical form. The results agree with experiment to within ~ 5%, 
for Eq = 300Mey, mg Kq — ISOMeF, and the hadronic scale parameter 
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(5 in the range 1.5 < /5 < 2.0. 

Sect. 4 outlines the formulation of the temperature dependent DSE (T- 
DSE for short) within the real time formalism [22] , instead of the imaginary 
time formalism a la Matsubara [23]. The order parameter in this regard 
may be chosen in one or more ways, a convenient choice being So which 
now "runs" with the temperature and is renamed as m^. Other analogous 
quantities which are expected to "run" with the temperature are the gluon 
mass renamed as m^t, and perhaps also the IR parameter whose connection 
with iTLt and irigt is brought out in Sect. (4). It is found that the both the 
constituent and gluon masses have the same "critical index" 7 = 1/3 (in 
accordance with the concept of 'universality' of critical indices), while the 
critical temperature works out at Tc 180 ± 2{)MeV. Sect. 5 concludes with 
a discussion including a comparison with contemporary approaches. 



2 Dyson- Schwinger Eq In Differential Form 

We start by writing the DSE in the Landau gauge which ensures that the A 
parameter does not suffer renormalization [24] . This is an additional precau- 
tion over and above the Pagels-Stokar DPT approach [10] to satisfy the WT 
identity. The starting DSE in the Landau gauge for the function S(p), after 
tracing out the Dirac matrices takes the form 



k^ - 2p.k 

1 2/x2 



^2 _|_ j^2 ^/j.2 _|_ 177,2^2 



The first term on the RHS corresponds to the rainbow approximation [9], 
while the second term gives the simplest realization of a gauge invariant 
structure by satisfying the WT identity a la Pagels-Stokar [10]. An analogous 
but slightly more involved ansatz due to Ball-Chiu [25] also can be seen to 
conform to the Landau form [10], through a visual inspection of both. To 
see this more explicitly, we list both forms for the relevant vertex functions. 
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first [25] (as given in [24]) followed by [10]: 



r.(p',p) = -t^A^ + A']/2 + f"^ _ \ [-^7.(p + p'){p, + p',)] + ^ _ (p, + p'^) 

z[p p ) p p 

r.{p\p) = -^7.+ ^^^^^^(P.+p'.) (2.2) 

where the momentum dependence {p,p') of the Ball-Chiu functions A, B is 
indicated by the unprimed and primed notations respectively and the mass 
function E(p) = B/A, while the Landau gauge corresponds to A = 1. The 
Ball-Chiu form [25] is seen to be compatible with Pagels-Stokar [10] (which 
is already in the Landau gauge A — 1). So, ithout further ado, we shall use 
only [10] for simplicity. 

We now adopt a suhtr active form of Renormalization by writing a similar 
equation for, say, p\ and subtracting one from the other. If p' is infinitely 
close to p, this results in a differential form. Thus we subject both sides of 
the eq.(2.1) to the differential operator p.d, not the scalar form p^dp2, since 
the former is more naturally attuned to handling two vectors p, k that occur 
on the RHS. The main advantage of this crucial step is to reduce the degree 
of divergence of the integral w.r.t. k by two units, which in turn allows 
further simplifications on S(p — k) on the RHS, since it falls off rapidly with 
k"^. In particular, we are allowed the following simplification as a result of 
this crucial step of reducing divergence via differentiation: 

F - 2p.k "^^'^^^^ 

A second simplification arises from a contraction of the factors {p — k)n{2p — 
k)i, and {Smuu — k^ki./k'^) which is almost independent of k^, and gives on 
angular integration [26] : 

2[p^ - {p.kf/k^] « 2(1 - n-^)p^; n = 4 

Further, against the background of the differential operator p.dp on both 
sides of (2.1), we can replace the mass function E^(p — k) inside the fermion 
propagator on the RHS due to an improved /^-convergence, by simply replac- 
ing Ij^(p — k) with S^(p), since this quantity already falls off with momentum 
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(see also [10]). The resulting eq.(2.1) now takes the form 



2E'(p) 




(p) - j:"{p) s(p) - s-(p)/2 

D{p-k) D^{p-k) 



(2.3) 



(4E(p)E'(p) + 2p2-2p.A;)][ 



k'^ + rrig {kP' + m^)^ 



where we have taken F1.F2 — —4/3, and defined derivatives and propagators 
as 



Note that decoupling of S(p) from k^ now facilitates the /^-integration, thus 
converting the DSE into a difi:erential equation, while the form of E(p) is as 
yet undetermined. (This structure is different from a more conventional one 
for a differential form of the DSE, by making the D(p — k) separable in terms 
ofp> andp<, etc [7, 14]). 

The next task is to integrate w.r.t. d'^k which for the o.g.e. term is 
still logarithmically divergent and hence requires 'dimensional regularization 
' (DR) a la t'Hooft-Veltmann [27], while the IR term gives a convergent in- 
tegral. We hasten to add that this divergence (despite the Landau gauge) 
may well be an artefact of the approximation E(p — /c) ~ S(p) in the nu- 
merator on the RHS, but since the divergence is only logarithmic, it is not 
sensitive to DR [27], and in any case it is a small price to pay for the huge 
advantage accruing from this (new) differentiation method for renormaliza- 
tion. Another approximation concerns the factor g"^ on the RHS of eqs(2.1) 
and (2.2) which, strictly speaking, is a function of the momenta p, k, but at 
this stage we must "freeze" the value of ag at a fixed value (to be specified 
below) so as to get a self-consistent asymptotic solution in the p^ = 00 limit. 
[More general solutions with the differential form (2.2) and variable ag have 
not been attempted here]. 

2.1 Dimensional Regularization for Integrals 

Denote the two integrals of Eq. (2.2) containing the o.g.e. term only by / and 
// respectively, of which only I is divergent (see above) , but // is convergent 



S'(P) 
D{p - k) 



(l/2)p.9pS(p)=p- 
T.\p) + {p-kf 



'%2^py, 



(2.4) 
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by itself. .Thus write for / in the Euchdean notation for dimension n, using 
the DR method [27, 11] 



, = 4,K2E'W - E»)C< / p,)„^(^_,,)(,..^„.) (2.5) 



rl. rco 

= ^gl{2Y.'{p)-Y."{p))C du dk'k 

Jo Jo 



2 un-2 ^ 



r(n/2)(27r)'*(A„ + fc2)2 



where we have introduced the Feynman variable 0<M<l,Cisa UV 
dimensional constant, e = 4 — n, and 

A„ = uE'^ip) +p'^u{l -u) + ml{l - u) 

The integration over A;^ is now straightforward , while that over u is simplified 
by dropping the rrig term since there is no infrared divergence.The result of 
all these steps after subtracting the UV divergence [27] is (with — Anag): 

I = (a,/7r)(2S'(p) - S"(p))[ln47r - 7 + 1 + HCVA,)] (2.6) 

where 

A,^^\p)+py2 (2.7) 

The other integral // which is UV convergent, does not need DR [27] and 
gives 

//^-K/.)(E(p)-E>)/2)(13Em±i^!) (2.8) 

Zip 

Thus the resulting DSE may be expressed compactly from (2.3) as 

2E'{p) = I + II + I' + II' (2.9) 

where we have taken / and II from the o.g.c. contributions (2.7) and (2.8) 
respectively, as well as added two similar (infrared)) terms /' and II' aris- 
ing from the second (IR) part of the gluon propagator (1.1). In the same 
normalization as above the last two work out as 

I' = ^(2E'(p)-E"(p))[ln(V^;)-l] (2.10) 
II' = -^(E(p)-E'(p)/2)(4E(p)E'(p)+/)[ln(V^J)-2] 



where we have made use of the smallness of TUg is simphfying some of the 
integrals over u. Note that the last two terms are at least of two lower orders 
in p than their o.g.e. counterparts, so that they will not contribute to the 
— oo limit of the differential equation (2.9). 

2.2 Large and Small Limits of DSE for T = 

To solve eq.(2.9), we try the ansatz [1, 20]: 

E(p) = Eo/[l + /5p2] (2.11) 

whose asymptotic form is compatible with pcrturbativc QCD expectations 
for massless quarks in the chiral limit [10]. And take the fixed value of p'^ in 
the argument of as at p^ = where is the UV parameter corresponding 
to the upper hmit of allowed in the solution of the DSE. [Other options 
exist, but not particularly convenient]. 

2.2.1 Lcirge limit 

Remembering the definition (2.4) for E', etc., we have in the large p^ limit 
for the function (2.11), 

E(p) ^ -S'(p) ^ +E"(p) 

Now remembering the upper limit of p'^ being constrained by the UV param- 
eter C^, substitution from (2.9) yields the result 

-2 = a,/7r[-3(ln(47r)-7 + l + ln2)-3]; 7r/a, = ^ MC/Ag) (2.12) 

Aq — 150MeV being the usual QCD scale parameter. Thus eq.(2.12) de- 
termines the value of the maximum momentum ( within this approach, and 
shows that our formalism does not permit p"^ to exceed C^. Unfortunately 
eq.(2.12), which corresponds to the check point p^ = oo, restricts C to a 
rather low value : 

C/Aq = 1.5490; C = 0.706GeVonly (2.13) 
where the MS scheme (not MS [28] ) has been employed. 
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2.2.2 Small limit 

Next we consider the small limit of eq.(2.9) where for the (fixed) argument 
of as we continue (for mathematical consistency) to maintain the same value 
of as corresponding to p"^ = C"^, leading after straightforward simplifications 
to 

Co + lnxi/xo-3+l/xo + [I' + II'] = 91n(C/AQ); Cq = hi 4% - = 1.9538; 

(2.14) 

where the dimensionless quantities are defined as 

xi = CP; xo = El(3 (2.15) 

Note that eq.(2.14) has a big term on the RHS, viz., 9 x 1.5490, needing 
a corresponding augmenting of the LHS, which can come only from the IR 
terms from (2.9), symbolically denoted by [/' + //'] in (2.14), that include the 
(as yet free) parameter A = 2/5/x^.[Of course these IR terms do not contribute 
to (2.12)]. 

2.3 Dynamical Mass And Mass Renormalization 

The third point p^ = —ml which defines the dynamical mass, corresponds to 
the "pole' of the propagator Sf{p), so that 

E^(imo) = +ml > E'^{0) (2.16) 

It may be recalled that a distinction between the dynamical and constituent 
masses already exists in the literature. Thus in the notation of Domb [29], 
(pp 322-324), p and po correspond to mo and Eq respectively]. 
Substituting from (2.11) gives a cubic equation in ml: 

J:l ^ mlil - f3mir (2.17) 

which implies that Eq < mo. Using the dimensionless variables Xq = /3Eq 
and yo = (3ml, ^^^^ reduces to the cubic yo(l — yo)^ = Xq which has at least 
one real solution for yo in terms of xq: 

13ml = yo^ 2/3 + ^[xo/2 - 1/27 ± ^Jxl/4 - Xo/27Y/' (2.18) 
± 
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whose nature can be seen as follows. For small xq, yo is also small (seen 
directly from the cubic form), but as Xq increases, yo increases more rapidly, 
until xo reaches a critical value Xc — 4/27 (seen from eq.(2.18)). Beyond this 
point Uo increases more slowly with xq. The corresponding critical value of 
p is 

/5e = 4/(27S2 ^ lM6GeV-\ /orSo = SOOMeV (2.19) 

We shall keep Sq fixed at SOOMeV, but vary j3 in the typical hadronic range 
1.0 < P < 2.0 for applications to key QCD parameters like < qq > and 
Now the propagator may be written as 

= Z... s'4"^;-f (2.20) 

making use of cq.(2.17), and formally introducing a "mass renormalization" 
factor Zm to be determined. However, using the condition (2.16) in the 
numerator and denominator of (2.19) shows immediately that near the pole 
the RHS already has the correct structure (mo — i'y-p)/[ml + p^], which 
suggests that = 1 ! On the other hand an alternative way to extract 
the factor (p^ + ml) from (2.19), suggests a non-zero value of Z^- This is 
seen from rewriting the RHS of (2.20) as 

Zmi^ip) - i-i.p)/[i:\p) - Y?{im) + + 

and taking the limit — > ^"^^o after extracting the factor (jn^ + p^) from 
the denominator. Z^ is now determined by the condition that at the pole 
this quantity reduce exactly to l/(mo + i^-p)- This gives 

Z^=(l-3m2/3)/(l-m2/3) (2.21) 

In view of this ambiguity in the working definition of Z^, it is not clear if this 
(finite) Zm is significant beyond unity. However within this subtractive renor- 
malization approach to the DSE, the divergences are already toned down to 
the logarithmic level, so that renormalization is probably less significant than 
for the usual (unsubtracted) DSE form. For the rest therefore we shall set 
Zrr, = 1 in what follows. 
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2.4 Solution of Eq.(2.9), Including IR Terms 



Taking account of the IR terms in (2.9) , the full equation (2.14) reads: 

= Axl-xo{B\ + l) + CX (2.22) 
A = -ln47r + 7 + 3-21n(C/So) + 91n(C/AQ) 
B = 7-61n(Eo/m<,); C = 2 - 2 ln(Eo/mg) (2.23) 

where xq = /9Sg, and A = 2^^j3. A practical way is to solve (2.22) for A 
with Eq = 300Mey and nig = 150MeV. This gives A for typical values 
of the 'range parameter p. Note that the connection between xq and yo is 
already determined by (2.17-18). Now with E = SOOMeV, a typical value 
P = IGeV-^, {xo = 0.09 and yo = 0.115), yields A = -0.0640. The latter 
is an index of the strength of a (small) IR term needed to provide a self- 
consistent solution of the DSE in the low momentum regime to match its 
solution for 'large' momenta. We shall come back to these quantities in 
Sect .4 for the T-dependent DSE. 

At this stage it may be asked as to what happened to the third check- 
point = —rriQ for the DSE, analogously to the points = oo and = 
considered in the foregoing. As a matter of fact, this condition has already 
been subsumed in the determination of the relation between the constituent 
and dynamical masses in eqs.(2.17-2.18) within the specific structure (2.11), 
so no new results can be expected from the DSE for = —ml. The check- 
point p^ = —ml will however come into play again in Sect. 4, but in a T- 
dependent form of the DSE. But before implementing the T-dependent DSE 
programme, it is first necessary to carry out two vital tests of this T = 
formalism, viz., its performance on the two crucial quantities < qq > and 
[1,5, 6] which we consider next. 

3 Tests of Mass Function: < qq > And 

The quark condensate and the pion decay constant are regarded as fairly 
sensitive tests of the mass function E(p) determined as a solution of the 
DSE, expressed in the differential form (2.9). To that end we first collect 
their formal definitions as follows. The condensate after tracing out the 
Dirac matrices is 

AN, /■ 4 E(p) 
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which simphfies on making use of eq.(2.19-2.20) to 



< (27r)4 y ^ ^ (p2 + mi) [(1 + - 2yo(l - Z/o)] ' ^'^^ 

where 

x^p^p; yo = "^o/3; Xo = Eg/? (3.3) 
The corresponding quantity may be defined in the chiral hmit by 
TV f 

VlP, = ^ j d'pTrmp^) + E(p2))75^Fi^(Pl)^7/.75-5Fi^(-P2)] 

where S'ir/j is given by (2.19-20), and P = P1+P2, and the pion-quark vertex 
function has been taken as [10] [S(pi)+S(p2)]/2/7r]. Fortunately the complete 
expression may be taken over from ref. [10] (also given in [1]), viz., 

f'- fd%. [1 - (PV4)g,.]S^(p) 

- (2^ y ^ (E2(p)+p2)2 (3-4) 

in the Euclidean limit. The derivations of (3.2) and (3.4) are shown in Appen- 
dices A and B respectively. The final result for the condensate is summarized 
in (A.4-A.6) where the standard Table of integrals [30] has often been em- 
ployed. Similarly, for the pion decay constant, the final result is given by 
(B.5). 



3.1 Results for Condensate And Pion Decay 

The key parameters of this theory are Eq, the constituent mass, and /3, the 
parameter for the non- perturbative hadronic scale. A third quantity, the 
dynamical mass mo is determined by these via cq.(2.17), which can be ex- 
pressed in terms of the dimensionless parameters Xq) and y^. Since the object 
of this investigation is not to provide a detailed phenomenological fit to these 
quantities, rather to see if this new differential form of the DSE is consis- 
tent with the conventional range of values of the constituent mass, we shall 
refrain from any fine-tuning and offer some typical values within this alter- 
native DSE framework, which is constrained by the fairly rigid connection 
between Eg and mo brought about by the cubic equation (2.17). Thus, with 
a fixed Eq at 300Mey, Table 1 depicts some typical values of {3, Xq and Uq 
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Table I: Variations of xq, yo With f3 





Xq 


yo 


1.00 
1.646 
2.00 


0.090 

4/27 

0.135 


0.115 
4/3 
1.365 



For these 3 sets we get under the MS scheme [28] 

< qq >o= (0.1545; 0.0932; 0.114)Eo//? (3.5) 

where we have depicted the sensitivity of this quantity to the main parameters 
P for Eo fixed at 300Mey . For the values hsted in Table I, the numbers work 
out at 

(359MeVf; {279MeVf; {28AMeVf 

respectively, suggesting that (3 should lie fairly close to its 'critical' value 
(3c = 1.646, without further tuning. Similarly, the pionic constant works out 
for the 3 values of Xq given above, as 

/2 = {92.0MeV)^; {93.1MeV)'^; {94.3MeVf (3.6) 

respectively, with Eq = 300Mey. This quantity is not sensitive to P but 
varies as the square of Eq. These values give a rough test of this formal- 
ism without vastly extending the numerical framework. Note that the the 
IR parameter A at —0.064 has been rather passive in these determinations, 
but its temperature dependence is going to play a more active role in the 
T-dependent DSE, for a self-consistent determination of the critical temper- 
ature Tc to be considered in Sect. 4 to follow. 

4 T-DSE In Real Time Formalism 

As noted in Sect. (1.1), since our DSE formulation departs from the more 
conventional thermodynamic formulations [20, 15] based on the free energy 
[15] or effective potential [20], we are not yet in a position to offer a complete 
set of critical indices near Tc, except the one for the T-dependence of the order 
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parameter rrit. Keeping this in mind, to formulate the T-dependent DSE, we 
have two broad options: real [22] vs imaginary [23] time formalisms. The 
T = structure of the DSE suggests that it is natural and convenient to 
employ the real time formalism and follow the prescription of Dolcn-Jackiw 
[22] for adding to the quark and gluon propagators (which can be easily 
read off from the main DSE, eq.(2.9)), the T - dependent imaginary parts of 
the Bose / Fermi types, leading to the modified propagators respectively as 
follows 



— ?■ 27r I ~ 

DFT{k) = ——^ + ^ 5{k' + ml); uj^Jml + k^ (4.1) 

k^ + m^ expu/T—1 ^ \ " 

Sft{p) = ^ ^ . , {p)+P) 4.2 

E(p) + I'y.p exp Ep/T + 1 

where the quark energy Ep is the fermionic analog of the gluon energy u, 
cq.(4.1). Taking the gluon case first, there are now two kinds of operations on 
(2.9). Namely, since the values arc being considered on the mass shell, we 
shall now write = — mf (instead of — mo)to emphasize the T-dependence 
of this quantity. Similarly (see Sect.l) wc shall consider the gluon mass mg 
and the constituent mass Sq to "run" with T, and designate them as Ht and 
mgt respectively. Considering the Bosonic and Fermionic Boltzmann factors 
(4.1-4.2)in this order, we shall have extra contributions to the four pieces on 
the RHS of (2.9), but giving rise to 3D integrals only. We now collect these 
values separately, first upgrading the T — results of Sect 2 to T 7^ 0. 

4.1 T-Dependent /; // and I';!!' 

To simplify the 4 pieces of the DSE, eq.(2.9), on the T-dependent mass shell, 
the following results are useful. 

S(p)=m,; 2S'(p) - S"(p) = (4.3) 

E(p) - E'(p)/2 = mt-(^ (4.4) 

Collecting these results on the (now T-dependent) 4 pieces on the RHS of 
(2.9) we have 

I + I' ^ ±M [2.9538 + In (2(7^^ + ^ [In {m'j2m%) - 1]I4.5) 

Ej mffJ 
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// + //' = (,„,-^(2-^)[l + -^[lnKV2m^J-2]] 

To these pieces must be added the T- parts of the gluon propagators 
(Bosonic) accruing from (4.1), and the T parts of the quark propagators 
(Fermionic) from (4.2). These are basically 3D integrals because of the 5- 
functions. To evaluate them the following quantities come into play 

D{p -k) = S=^(p) + (p - kf = -mlt + 2mtUJ (4.6) 

4E(p)E'(p) + 2p2 - 2p.k = - 2ml + 2^^^; (4.7) 

Et 

Here we have taken the rest frame of p^, viz., p = 0. The Bosonic T- parts 
normalized to the pieces in (4.5) are 

BOSEt = 4 / - ml [Ibt + Hbt] (4.8) 

J ^ " expcu/T — 1 

where the lower limit of cu integration is nigt and the two integrands are 



(4.9) 



Similarly for the Fermionic parts, denoted by I(FT) and II(^FT) respectively. 
The complete T-dependent DSE is then obtained by modifying (2.9) a la (4.5) 
and adding the pieces (4.9-4.10), and the corresponding Fermionic parts, after 
integrations. Before carrying out these integrations we notice some general 
features of these quantities in the neighbourhood of the critical temperature 
Tc. Namely, 

i) the powers of rrit are spaced by three units; 

ii) rrit and rrigt are always involved in identical ratios. 

One may infer from this that the critical index 7 for both is the same at 
37 = 1, consistent with universality [29] for such quantities. Thus in the 
neighbourhood of Tc one may take 

K; mgt] ^ [So; m,]T^; r = 1 - T/T,; 7 = 1/3 (4.11) 
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The the T-dependence of mo may be handled via (2.17). Note that in the 
neighbourhood of Tc, Ri rrit, a result which is consistent with eq. (7.278) 
of ref [29]. Retaining only the lowest powers of the small quantities m^, m^t, 
most of the terms in the T — DSE will drop out, and the integrals over 
(4.9-10) will lead to the net bosonic contribution 

To this T-dependent (gluon propagator) contribution, must be added the 
corresponding quark propagator contribution, eq.(4.2), near T = Tc, by fol- 
lowing a similar procedure to above. For brevity, we indicate only the extra 
features, before writing the final result. The Fermionic T- part of the quark 
propagator in(4.2) now becomes 

{exp{E{p-k/T) + l)) 

And analogous to (4.6), 



k"^ + ml ^ 2mtEk - 2m'^t + mlt, E^ = + m^) (4.14) 

Next, taking account of eqs. (4.3-4), and proceeding as in the gluon case, 
we can evaluate the quark counterpart of eq.(4.8) in the neighbourhood of 
T = in the form 

FERMIt « [-r/?m,2tan-^[— ^— ] + At(-7 + ]n{T/mt))/i]/{in^) (4.15) 

1 +mt 

It is easily checked that this quark contribution is at least of 0{y/j3mt) com- 
pared with the gluon's, so that it is justified to neglect it, at least near the 
Critical point. The master equation T — DSE, keeping only the lowest order 
terms, now simplifies to 

AX, 

Li + BOSEt = 0; Li = In (m^/2m5) - 2 (4.16) 



(3mt 



This equation suggests a simple structure for Aj, perhaps one of the few that 
are consistent with its solution, viz., 

A* = Xo{mgt/mgy[-\nT^ + 1]; r = 1 - T/T, (4.17) 
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where Aq may be identified with the value found in Sect. 3, viz., A = — 0.064± 
0.003, and the term unity in square brackets signifies its normahzation at T = 
0. Eq.(4.16) after substitution from (4.12), now reduces to two equations, 
involving the coefficients of 

t'^] r'^lnr''' 

respectively, but we skip these equations for brevity. The result, after elimi- 
nation of the quantity Li of (4.16) from them, and dividing out by T, is 

-1/2 + (1/21/) [1 - 1/72] = l/21n[27;/m^]; v = m^/Eo (4.18) 

using (4.11) near the Critical Point. Substituting from Sect. 3. 3, viz., i/ ^1/2 
gives the surprisingly simple result 

i/ftil/2; 2Tc^m3Cxp7/8 (4.19) 

leading to a reasonable value for the critical temperature, viz., 

Tc ^ 180 ± 2QMeV] {nig = 15QMeV). (4.20) 

4.2 Condensate and pionic constant near T = Tc 

For completeness we offer some brief comments on the predictions of this 
simple formalism on the corresponding T-dependent quantities < qq > and 
near the critical point T = Tc, analogously to the results of ref [1,5]. This 
is possible in view of the analytical expressions for these quantities as given 
in Appendices A and B in terms of yQ,a and Xq respectively. In T-dependent 
form, yo ~ mf,a ~ m^, and Eq ~ ''^t, which in turn are expressible in 
terms of the basic 'order' parameters mt and nigt, eq.(4.11). Substitution in 
eqs.(A.6) and (B.5) shows that < qq > and tend to zero neat the Critical 
Point like and Inl/nif respectively, in general accord with standard 
expectations. 

For comparison with other approaches, the chiral perturbation theory [6] 
for /^T) predicts [5] 

Nf r 

MT) = /.[I - / dMew (E/T) - l]-'/E] 

which however does not indicate how this quantity behaves near T^. Another 
form [1] which is more in line with our parametrization of S(p), suggests that 
Et should vary as < gg >t, in agreement with our result for the condensate. 
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5 Summary and Conclusion 



In retrospect, we have proposed a new (differential) form of the Dyson- 
Schwinger Equation (DSE) for the mass function based on an (in- 

finitesimal) subtr active form of Renormalization in QCD. Such 'subtraction' 
in turn amounts to employing a differential operator of the form p^(9^ apphed 
on both sides of the DSE, whose effect on the RHS is to lower the degree 
of divergence w.r.t. the integration variable A-^j by TWO units. It is in the 
background of this (differential form of) subtractive renormalization, that it 
becomes possible to approximate the quantity S(p — /c) inside the integral by 
S(]?) since the A;-dependence of this already decreasing quantity is no longer 
sensitive to the principal term {p — kY in the quark propagator. [Without 
this background of an improved A;- convergence, however, this approximation 
would not have been justified ]. This crucial step which has facilitated the 
integration over d*k without further ado, has thus helped convert the DSE 
into a second — order differential equation, the extra order (beyond the rain- 
bow approximation [9]) arising from the term responsible for satisfying the 
WT identity a la Pagels-Stokar [10], so as to preserve gauge invariance. To 
reinforce this effect, we have employed the Landau gauge which makes the 
DSE virtually depend only on the mass function S(p) by effectively elimi- 
nating the A-function [27]. (The 'ghost terms' do not appear in this effective 
description) . 

To solve the resulting differential form of the DSE, we have taken recourse 
to three crucial check-points: = oo, = 0, and = — mg, using a pole 
ansatz, (2.16), (c.f. [1]) which is consistent with the form p~^ in the large 
p2 regime, in agreement with dynamical breaking of chiral symmetry for 
massless quarks [10], provided the argument of as is held fixed at some chosen 
value (here the UV parameter Q . This has given a rather small value on the 
UV parameter ( that appears as an argument of as, which effectively restricts 
the range of applicability of this formalism to moderate values of p"^ (perhaps 
adequate for the cosmological application envisaged in this paper). For the 
low p^ regime, we have introduced two kinds of masses: the constituent mass 
E(0) which is generally believed to be of ~ SOOMeV, and the dynamical 
mass mo which satisfies the equation E(imo) = tuq corresponding to the 
pole position of the quark propagator = —ml (see also [29]). Now for 
the simple form (2.16) the connection between the two 'masses' is given by 
(2.18) (as the solution of a cubic equation) which corresponds to E(0) < 
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mo. The parameter (3 in eq.(2.16), for a given S(0), has been taken as a 
typical hadronic scale befitting the low energy regime of the DSE. The small 
IR parameter 2/x^(= A//3) which has played a passive role in the T = 
description of the DSE, turns out to be rather crucial for T > 0, for which 
the ansatz (4.17) is necessary for a self-consistent solution of the T — DSE 
(see further below). We have also considered a non-zero value of the gluon 
mass for which several arguments have been advanced in the literature [13- 
15]. 

We have also carried out two important applications of E(j9) obtained 
from this new formulation of the DSE, viz., the quark condensate and the 
pion decay constant, more by way of some basic calibration of the formalism 
than as a means of detailed phenomenological fits to hadronic data. Thus a 
fit to within < 10% has helped fix the parameters involved. After this check, 
we have attempted in Sect. 4 a T-dependent formulation of the DSE to see 
the extent to which it can simulate the critical temperature and at least one 
of the critical indices. To that end, we have taken the — —ml limit of 
the T-DSE near the Critical Point Tc where it is small. In this respect, the 
demands of consistency have necessitated a T-dependence of the IR confining 
parameter A, for which an ansatz of the form (4.17), calibrated to its value 
at T = 0, is indicated. Two clear results have emerged from the analysis, 
viz., i) a bunching of the powers of in units of three suggest a critical 
index 7 = 1/3 according to conventional analysis [29] ; ii) and the 'matching' 
of the coefficients of hke powers of the reduced temperature r have led to 
a very simple solution of the form (4.18), leading to the reasonable Tc at 
180 ± 20Me\/. 

For a comparison of this result with those of contemporary approaches 
[20,15], our approach differs from these in one important respect: the role 
of an external if-field is sought to be partially simulated by the IR param- 
eter A which is necessarily T-dependent, instead of by small but non-zero 
u — d masses [20,15]. Further, in view of our explicit analytical expressions 
for < qq > and we have also obtained analytic structures for their T- 
dcpcndence, and found indeed that they both vanish at the critical point, 
without a detailed numerical analysis [6,15, 20]. However this approach has 
its weak points, especially the ad hoc nature of eq.(4.17) for the T-dcpcndence 
of the IR term. A second one is lack of a more plausible understanding of 
the extent to which the IR term can substitute for the current masses [20, 
15] to simulate the if-field effect. Attempts at throwing more light on these 
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issues, as well as extending the T-DSE formalism to facilitate the evaluation 
of other critical indices [20, 29], are envisaged. And in view of its central 
role, several other applications of the "mass function", such as tt — > 27, and 
e.m. pion form factor at finite temperature [31], are under way. 
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7 Appendix A: Evaluation of < qq > 







Using the notations x = and ijq = Pml, and anticipating a UV divergence 
which requires a DR treatment [27], we write 4 ^ n in eq.(3.2) which reduces 
after the angular integration [27, 11] to 



(27r)«r(n/2)/3"/2-i Jo 
where 

^i"^) - r^4rr£K-2 — st' = 2t/o(i - yo) (a.2) 

{x + yo)[{l + xy - 
Now break up F{x) into partial fractions 

1 A-yo l-yo + a 1 -yo - a . 



F{x) = 



[{1-yoy -a^]'x + yo 2(1 + a) 2{l + x + a)' 
The integration of each term above is carried out according to 

x'^/^-^dx/iA + x)^ A^-'/^r{n/2)r{l - n/2) (A.3) 



Jo 



The rest is a matter of collecting all the 3 terms after giving a DR [27] 
treatment to each. The final result is 

<qq>o = -^[^(a)(ln47r-7 + l + lnC'/3) + M«)]; 
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^ ^ (1 - yo^ - a2 ^ / ^ ^ 1 - yo ± a 

This result is valid for small yo (i.e., /3 ~ 1), when > 0. However for larger 
/3, vide eqs(2. 18-19) of text, yo exceeds unity, and < 0. For such cases, put 
— —b^. In particular the partial fraction break-up for /3c, corresponding 
to yo = 4/3, is rather simple: 

F{x) ^ + (1 + 1/3 



(1 + x)2 + 62 a; + 
since 6^ + (yo — 1)^ becomes unity. Now using the result [30] 

(1 + 62)n/4-i sin (^/2 - l)t 



f 

Jo 



10 l + b'^ + 2x + x'^ sintsinn7r/2 
where 



(A.5) 



cos t = +1 /Vl + sin i = +6/ VTTF 
and giving a DR treatment [27] as above the corresponding result to (A.5) is 



<qq>o = -^[6^(ln47r-7 + l + lnC'/?)+/(6)]; 

f{h) = (l/6-672)ln(l + 62) _lnyo/3- (4V3)tan-^6 (A.6) 



For purposes of obtaining the tcmpcraurc dependence of the quark con- 
densate (to be discussed in Sect. 4 of the text), we record the results of these 
integrations in the limit of small a and yo, for which eq.(A.4) is appropriate: 

51(0) ~ 1 -|- 0(yo); /i(a) ~ yolnyo -|- ~ alna (A. 7) 

Substitution in (A.l) gives in this limit 

< qq >o= [(In 47r - 7 + 1 + In CP) + 0(a In a)] (A.8) 

which lends itself immediately to a finite T treatment in the neighbourhood 
of the Critical point (see text, Sect.4). 
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8 Appendix B: Evaluation of 

Since the integral (3.4) is convergent by itself, DR [27] is not needed in this 
case. After the angular integrations (using the dimensionless units x, i/q as 
before ), and carrying out the differentiations, (3.4) reduces to 

where the integral is defined by 

;i + x)(i + 3x/2 ) 

[xo + x{l + xy] 
Now transform the variable from x to u, as 



^ = /„ dxx )^ , '^^^ , (B.2) 



u^^^; 0<u<l. 
The result of this is to give an integral in u as 

r^ duu{l-u){l + u/2) 

Jo [xo{l -u)^ + uY ^ ■ ^ 

While this integral is i principle exactly doable, it is instructive to obtain an 
approximate analytical expression which in practice is sufficiently accurate, 
so as to lend itself to a generalization to finite temperatures (see below). 
The trick lies is the observation that most of the contributions arise from the 
region of small values of u. Then (B.3) simplifies to 

duu{l — u/2) 



10 [xo{l — "iu) -\- uY 
Now integration by parts gives the final result 

Unlike the case of < >, this result is valid for all allowed Xq. For purposes 
of determining the temperature dependence of to be discussed in Sect. 4, 
we record as in Appendix A, the corresponding results for small xq This gives 

Tfa l/21nl/a;o - 1/2 
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Substitution in (B.l) leads to 



a ^iNc ri 1 / 1-3x0, 



which lends itself immediately to a finite T treatment in the vicinity of he 
Critical Point (see text, Sect .4). 
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